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Abstract 

In this paper, we prove the existence of non-negative solutions for a 
non-local higher order degenerate parabolic equation arising in the mod- 
eling of hydraulic fractures. The equation is similar to the well-known 
thin film equation, but the Laplace operator is replaced by a Dirichlet-to- 
Neumann operator, corresponding to the square root of the Laplace op- 
erator on a bounded domain with Neumann boundary conditions (which 
can also be defined using the periodic Hilbert transform). In our study, 
we have to deal with the usual difficulty associated to higher order equa- 
tions (e.g. lack of maximum principle). However, there are important 
differences with, for instance, the thin film equation: First, our equation 
is nonlocal; Also the natural energy estimate is not as good as in the case 
of the thin film equation, and does not yields, for instance, boundedness 
and continuity of the solutions (our case is critical in dimension 1 in that 
respect). 
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1 Introduction 

This paper is devoted to the following problem: 

lit + [u n I(u) x ) x = for x G fi, t > 
< u x = 0, u n I{u) x = for x edil, t > (1) 
u(x, 0) = uq(x) for x € £1 



"CEREMADE, Universite Paris-Dauphinc, UMR CNRS 7534, place de Lattre de Tassigny, 
75775 Paris cedex 16, France 

t Department of Mathematics. Mathematics Building. University of Maryland. College 
Park, MD 20742-4015, USA 



1 



where fl is a bounded interval in R, n is a positive real number and / is a non- 
local elliptic operator of order 1 satisfying I o I — —A; the operator / will be 
defined precisely in Section [3J as the square root of the Laplace operator with 
Neumann boundary conditions. When = R, it reduces to / = —(—A) 1 / 2 . In 
the sequel, we will always take ft = (0, 1). 

When n = 3, this equation arises in the modeling of hydraulic fractures. In 
that case, u represents the opening of a rock fracture which is propagated in 
an elastic material due to the pressure exerted by a viscous fluid which fills the 
fracture (see Section [5] for details). Such fractures occur naturally, for instance 
in volcanic dikes where magma causes fracture propagation below the surface of 
the earth, or can be deliberately propagated in oil or gas reservoirs to increase 
production. There is a significant amount of work involving the mathematical 
modeling of hydraulic fractures, which is beyond the scope of this article. The 
model that we consider in our paper, which corresponds to very simple fracture 
geometry, was developed independently by Geertsma and De Klerk and 
Khristianovic and Zheltov [33] • Spence and Sharp [JTJ initiated the work on 
self-similar solutions and asymptotic analyses of the behavior of the solutions of 
|T]) near the tip of the fracture (i.e. the boundary of the support of u). There 
is now an abundant literature that has extended this formal analysis to various 
regimes (see for instance [T], [2J, [33J and reference therein). Several numerical 
methods have also been developed for this model (see in particular Peirce et 
al. [3S], [37], [31] and [35]). However, to our knowledge, there are no rigorous 
existence results for general initial data. This paper is thus a first step toward 
a rigorous analysis of ([1]). 

From a mathematical point of view, the equation under consideration: 

u t + {u n I{u) x ) x = (2) 

is a non-local parabolic degenerate equation of order 3. It is closely related to 
the thin film equation which corresponds to the case / = d xx : 

u t + (u n u xxx ) x = (3) 

(note that the porous media equation corresponds to the case I(u) = u). In 
particular, like the thin-film equation, Equation ([2]) lacks a comparison principle, 
and the existence of a non-negative solution (for non-negative initial data) is 
thus non-trivial (it is well known that non-negative initial data may generate 
changing sign solutions of the fourth order equation dth + d xxxx h = 0). 

However, compared with ([3J, the analysis of presents some additional 
difficulties: First, the operator / is non-local and the algebra is not as simple as 
with the Laplace operator. Second, because of the lower order of the operator 
I, the natural regularity given by the energy inequality (u € H ? rather than 
u £ H 1 ) does not give the boundedness and continuity of weak solutions even 
in dimension 1. 

A remarkable feature of ([2]) and (J3j is that the degeneracy of the diffusion 
coefficient permits the existence of non-negative solutions. In the case of the 



2 



thin film equation ([3]), the existence of non- negative weak solutions was first 
addressed by F. Bernis and A. Friedman [TT| for n > 1. Further results were 
later obtained, by similar technics, in particular by E. Beretta, M. Bertsch and 
R. Dal Passo [5] and A. Bertozzi and M. Pugh [T31 Q3]. Results in higher 
dimension were obtained more recently (see [2"5I |2"TI l2"0] ). 

As in the case of the thin film equation, our approach to prove the existence 
of solutions for ([1]) relies on a regularization-stability argument, and the main 
tools are integral inequalities which we present now. 



Integral inequalities. Besides the conservation of mass, the solutions of © 
satisfy two important inequalities (that have a counterpart for the thin film 
equation): Assuming that we have fl = R and / = —(—A) 1 / 2 (for the time 
being), we can indeed easily show that the solutions of ([2]) satisfy the energy 
inequality 

i-T 



u n (I{u) x ) 2 dxdt < - / u Q I(u Q )dx (4) 



u(t)I(u(t))dx 



(where — J ul(u) dx is the homogeneous H 1 / 2 norm) and an entropy like in- 
equality 



G(u(t))dx 



u x I{u) x dxdt= I G(uo)dx 



(5) 



where G"(s) = ±. 

The energy inequality © controls the L°°(0, T; i/ 1 / 2 (fi)) norm of the so- 
lutions (instead of L°°(0, T; iJ 1 (fi)) for the thin film equation). We see here 
that the order 1/2 for the operator / is critical in dimension 1 in the sense that 
we are just short of a L°°((0, T) x fi) estimate and continuity of the solutions. 
Because of that fact, many of the arguments used in the analysis of the thin 
film equation do not apply directly to our case. 

Next, we observe that as in the case of the thin film equation, the entropy 
inequality ([5]) provides some control on some negative power of u for n > 2. 
Indeed, we can take 



G(a) 



■ dtdr 



so that G is a nonnegative convex function satisfying G'(l) = and G(l) = 0. 
This yields: 



sins 



„2-n 



G(s) = 



(2-n)(n-l) n-l 2-7i 



In — h s 
s 



(n- 2){n - 1) s"- 2 n-l ?i-2 



when n = 1 
when 1 < n < 2 
when n = 2 
when 7i > 2. 



(6) 
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Note that G(0) = +00 when n > 2, while G is bounded in a neighborhood of 
for n G [1,2). This will be key in proving the non- negativity of the solution. 
The entropy equality also gives some control on the L 2 (0, T; H 3 / 2 (n)) norm of 
the solution which will be crucial in getting the necessary compactness in the 
construction of the solution. However, in order to make use of this inequality, we 
need J n G(uq) dx to be finite, which, when n > 2 prohibits compactly supported 
initial data. 

Besides those two inequalities, there are several other integral estimates that 
have proved extremely useful in the study of the thin film equation. The simplest 
one are local versions of d3| and (|S|). However, because of the nonlocal character 
of the operator J, it seems very delicate to establish similar inequalities for @. 

Another crucial estimate in the analysis of established by Bernis [10], is 
the following: 

{u^fdx <C u n u 2 xxx dx 



for n G (5,3). Such an inequality yields important estimates from the dissipa- 
tion in the energy inequality (despite the degeneracy of the diffusion coefficient). 
Again it is not clear what would play the role of this inequality in our situa- 
tion. The same remark holds for the so called a-entropy [IT] [5] [14] [15]: for 
a G (max(— 1, | — n), 2 — n), a 7^ 0, it can be proved that the solutions of the 
thin film equation ([3]) satisfy: 

^-y J u a+1 (-,T)dx + C £ J (Idxu^^f + ld^u^^l^dxdt 

< — !— [ < +1 dx. 

These last two inequalities are essential in establishing many qualitative 
properties of the solutions, such as finite speed expansion of the support and 
waiting time phenomenon. Though we expect such properties to hold for j2]) as 
well, it is not clear at this point how to deal with the non local character of /. 

Finally, let us comment on the power of the diffusion coefficient u n . Inter- 
estingly, the power n = 3, which is the physically relevant power in our model, 
is critical in many results for the thin film equation. In particular many exis- 
tence and regularity results (as well as waiting time results) are only valid for 
n G (0, 3). It is actually believed that for n > 3, (and it is proved for n > 4) the 
support of the solutions of ([3]) does not expand. It is not clear what would be 
the critical exponent for @, though numerical results suggest that for n = 3, 
the support of the solutions of ([2]) docs expand for all time. 

Main results. A weak formulation of ([2]) is given by 

// udtipdxdt+ // u n d x I(u) d x (p dx dt = — / UQip(0,-)dx 
JJQ JJq Jn 
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for all <p e T>(Q) where Q denotes f2 x (0, T). However, because of the degeneracy 
of the coefficient u n , it is difficult to give a meaning to the term u n d x I(u). 
We thus perform an additional integration by parts to get the following weak 
formulation of ©: 



u dfif dx dt — J J nu n d x u I(u) d x ip dx dt — j j u n I(u) d xx ip dx dt 

= - / u <p(0,-)dx (7) 
Jn 

for all ip e T>(Q) satisfying d x tp\gn = 0. 

We are going to prove the following existence theorem: 

Theorem 1. Assume n > 1. For any non-negative initial condition uq € 

G(uq) dx < oo, (8) 



/ 



/Q 

there exists a non-negative function u £ L°°(0, T, H% (O)) smc/i £/iai 

ueL 2 (0,T,#|(Q)) (9) 

which satisfies (J7J) /or a// 93 £ satisfying d x (p\gn = 0. 

Furthermore u satisfies, for almost every t £ (0, T) 

u(t,-)dx = / uodx, (10) 

where the function g <G L 2 (Q) satisfies g = d x (u% I(u)) — ^u^~d x u I(u) in 
£>'(Q), and 

I G(u(x,t))dx + \\u\\ 2 3 < / G(u (x))dx. (12) 
Jn i J (o,t;H](f))) Jo 

We recall that the function G : (0, 00) —> K+ is given by (|6]). The space 

3 

Hjy(Q) appearing in (|9]) will be defined precisely in Section [3] In particular, the 
following characterization will be given: 

where d(x) denotes the distance to 3f2. Condition ([9]) thus implies that u 
satisfies u x — on <9f2 in some weak sense. 

Note that at least formally, we have g = u%d x I(u) (though we do not have 
enough regularity on u to give a meaning to this product in general). Finally, 
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we point out that we have H^(fl) C W l p (fl) for all p < oo and so every terms 
in ([7]) makes sense. 

For n > 2, condition §5§ requires in particular that Supp(uo) = Q and 
inequality (fT2|) implies that this remains true for all positive time. This is a 
serious restriction since the case of compactly supported initial data is physically 
the most interesting (see Section [2]). We hope to be able to get rid of this 
assumption in a further work. 

For n > 3, we can actually show that condition © requires u(-,t) to be 
strictly positive for a.e. t € (0,T). In fact, we can prove: 

Theorem 2. When n > 3, there exists a set P C (0,T) such that |(0,T)\P| = 
and the solution u given by Theorem [7] satisfies 

u(-,t) G C a (tt) for allte P and for all a<l, 

and u(-,t) is strictly positive in tt. Finally, u solves 

d t u + d x J = inV'(n), 

where 

J(-,t) = u n d x I{u) 6 L\n) for all t 6 P. 



Organization of the article. The paper is organized as follows: In Section[5J 
we give a brief description of the mathematical modeling of hydraulic fracture 
which gives rise to equation ([5]) with n = 3. In Section [3J we introduce the 
functional analysis tools that will be needed to prove Theorem [TJ In particular, 
the non-local operator I(u) is defined, first using a spectral decomposition, then 
as a Dirichlet-to-Neuman map. An integral representation for /, using the 
periodic Hilbert transform is also given. Section |4] is devoted to the study of a 
regularized equation while the proof of Theorem [TJ is given in Sections [5] (for the 
case n > 2) and |6] (for the case n S [1,2)). Theorem [2] is proved in Section [71 

Acknowledgements. The authors would like to thank A. Pierce for bringing 
this model to their attention and for many very fruitful discussions during the 
preparation of this article. The first author was partially supported by the ANR- 
projects "EVOL" and "MICA" . The second author was partially supported by 
NSF Grant DMS-0901340. 

2 The physical model 

When n = 3, Equation @ can be used to model the propagation of an im- 
permeable KGD fracture driven by a viscous fluid in a uniform elastic medium 
under condition of plane strain. More precisely, denoting by (x,y,z) the stan- 
dard coordinates in R 3 , we consider a fracture which is invariant with respect 
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to one variable (say z) and symmetric with respect to another direction (say 
y). The fracture can then be entirely described by its opening u(x,t) in the y 
direction (see Figure [2]). Since it assumes that the fracture is an infinite strip 
whose cross-sections are in a state of plane strain, this model is only applicable 
to rectangular planar fracture with large aspect ratio. 

We now briefly describe the main steps of the derivation of ([T]). 



2.1 Conservation of mass and Poiseuille law 



The conservation of mass for the fluid inside the fracture, averaged with respect 
to y yields: 

d t (jm) + d x q = inK (13) 

where p is the density of the fluid (which is assumed to be constant) and q = 
q(x, t) denotes the fluid flux. This flux is given by 

q = puv (14) 

where v is the y-averaged horizontal component of the velocity of the fluid 



1 



u/2 



u/2 



VH{t,x,y) dy. 



Under the lubrification approximation, Navier-Stokes equations reduce to 

d 2 v H 



/'■ 



dy 1 



-(t,x,y) = d x p(x,t) 



where p denotes the pressure of the fluid at a point x and p denotes the viscosity 
coefficient. Assuming a no-slip boundary condition v — at y = ±u/2, we 
deduce 



v H (t,x,y) = -d x p 



1 2 1 2 



for — u < y < u 
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and so 



v(x,t) = -—d x p(x,t) 
12/i 



Using (|14p . we deduce Poiseuille law 



q = -T^ d * p - (15) 

Together with (fT3|) . this implies 

In order to obtain (|3J| , it only remains to express the pressure p as a function of 
the displacement u (i.e. p = —I(u)). 

2.2 The pressure law 

For a state of plane strain, the elasticity equation expresses the pressure as a 
function of the fracture opening. After a rather involved computation, one can 
derive the following nonlocal expression (see [T5]): 



4?r J r z-x 



where E' denotes Young's modulus. Denoting by T~L the Hilbert transform, we 
can rewrite this formula as 

P (x,t) = ^-H(d x u) = IVa^VM) 

where (—A) 1 / 2 is the half-Laplace operator, defined, for instance, using the 
Fourier transform by J r ((— A) 1 / 2 ?/) = |£|.F(u). 

It is well known that the half Laplace operator can also be defined as a 

Dirichlct-to-Ncumann map for the harmonic extension. More precisely, the 
pressure is given by 

p(x) = ^d y v(x,0) (16) 

where v solves 

-Av = in K x (0, +oo), , , 

v(x,0) = u(x,t), on R. ^ > 

By taking advantage of the symmetry of the problem, the function v(x,y) can 
be interpreted as the displacement of the rock. Denoting I(u) = — (— A) 1 / 2 ^, 
we deduce p = —^-I(u) and so 

dtU+^-d x (u 3 d x I(u)) = 0. 
4o^ 
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A technical assumption. In order to reduce the technicality of the analysis, 
we will assume that the crack is periodic with respect to x. Since we expect 
compactly supported initial data to give rise to compactly supported solutions 
whose supports expand with finite speed, this is a rcasonnablc physical assump- 
tion. We also assume that the initial crack is even with respect to x = and 
we look for solutions that are also even. 

By making such assumptions (periodicity and evenness), we can replace (jf 7|) 
with the following boundary value problem 

-Av = in x (0,oo) 

v v = on dfl x (0, oo) 

v{x, 0) = u(x) on O 

with il = (0, 1) if the period of the initial crack is 2. The cylinder f2 x (0, +oo) 
is denoted by C in the remaining of the paper. 

Mathematical definition of the pressure. It turns out that it is easier to 
define first the operator / by using the spectral decomposition of the Laplace 
operator: We take {Xk,fk} the eigenvalues and corresponding eigenvectors of 
the Laplacian operator in fi with Neumann boundary conditions on dfl: 

( -Aip k = X k ip k in tt 
\ d u (flk =0 on dfl. 

We then define the operator / by 

oo OO ^ 

I(u) := ^c k ip k {x) i ^ -y^CfcA^yfc(x) 

k=0 k=0 

which maps i? 1 (SI) onto L 2 (f2). We will prove that this operator can be char- 
acterized as a Dirichlet-to- Neumann map (see Proposition [2]) and that it has an 
integral representation as well (see Proposition [3} . 

2.3 Boundary conditions 

The opening u is solution of @ in its support. The model must thus be sup- 
plemented with some boundary conditions at the tip of the fracture. 
Assuming that Supp(it(i, •)) = [—£(t),£(t)}, it is usually assumed that 

u = 0, u 3 d x p = at x = ±£{t) 

which ensures zero width and zero fluid loss at the tip. 

We point out that the model is a free boundary problem, since the support 
[— £(t),£(t)] of the fracture is not known a priori. Since the equation is of order 
3, those two conditions are not enough to fully determine a solution. In fact, 
there should be an additional condition which takes into account the energy 
required to break the rock at the tip of the crack. Consistent with linear elastic 
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fracture propagation, we can assume that the rock toughness Kj c equals the 
stress intensity factor Kj. When the crack propagation is determined by the 
toughness of the rock, a formal asymptotic analysis of fracture profile at the tip 
(see [J [53]) then shows that 

X' 

u(x, t) ~ — s/W) - as x -> ^(*) ( 18 ) 

with K' = 4yJ^Kj c (and a similar condition for x — > —£). One can now take this 
condition on the profile of u at the tip as the missing free boundary condition. 
The resulting free boundary problem is clearly very delicate to study (remember 
that ([T]) is a third order degenerate non-local parabolic equation). 

A particular case which is simpler and still interesting is the case of zero 
toughness (Ki c = 0). This is relevant mainly if there is a pre- fracture (i.e. the 
rock is already cracked, even though u = outside the initial fracture). Math- 
ematically speaking, this means that Equation ([I} is now satisfied everywhere 
in R even though u is expected to have compact support. No free boundary 
conditions are necessary. This is the problem that we are considering in this 
paper. 

Note that one should then have lim x _^(i(t) — x)~ 1 / 2 u(x, t) = at the tip of 
the crack. In fact, formal arguments show that the asymptotic behavior of the 
fracture opening near the fracture tip should be proportional to (£(t) — a;) 2 / 3 
(see [U[23]). 

This approach is very similar to what is usually done with the porous media 
equation, and it has been used successfully in the case of the thin film equation 
to prove the existence of solutions with zero contact angle (in that case, we speak 
of precursor film, or pre- wetting) . The study of the free boundary problem with 
free boundary condition (|18[) would be considerably more difficult (one would 
expect the gradient flow approach developed by F. Otto [35] for the thin film 
equation with non zero contact angle to yield some result when n = 1). 



3 Preliminaries 

In this section, we define the operator / and give the functional analysis results 
that will play an important role in the proof of the main theorem. A very similar 
operator, with Dirichlct boundary conditions rather than Neumann boundary 
conditions, was studied by Cabre and Tan [17]. This section follows their anal- 
ysis very closely. 



3.1 Functional spaces 

The space H s N (Vt). We denote by {Afc, (pk}k=o,i.2... the eigenvalues and corre- 
sponding eigenfunctions of the Laplace operator in Q with Neumann boundary 
conditions on d£l: 

-Aip k = \ k tp k in fl , , 

d v Lpk. =0 on dfl, 
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normalized so that J^tp^dx = 1. When £7 = (0, 1), we have 

Ao = , ipo(x) = 1 

and 

Afc = (kir) 2 , tfk(z) = V2cos(knx) fc = 1, 2, 3, . . . 

The cpk's clearly form an orthonormal basis of L 2 (f2). Furthermore, the ipk's 
also form an orthogonal basis of the space iTjy(fi) defined by 



equipped with the norm 



fe=0 fc=0 



k=0 

or equivalently (noting that cq = IMIl 1 ^) and > 1 for k > 1): 



2 _ ll„,l|2 , II ,,||2 



r (n) - \m\LHn) + WIh^q) 



where the homogeneous norm is given by: 

oo 

IHIirj,(f2) = ^2 c lK- 

k=l 

A characterisation of i?jy(f2). The precise description of the space H^(il) 
is a classical problem. 

Intuitively, for s < 3/2, the boundary condition u v = does not make sense, 
and one can show that (see Agranovich and Amosov [4] and references therein) : 

F^r(O) = H s (fl) for all < s < |. 
In particular, we have H^(il) = H 2 (f2J and we will see later that 

2 / /„./„.\ „.f„.\\2. 



" U "ii3(fi) = J J ~ u ( x )) "(x,y)dxdy 

where y) is a given positive function; see ((23)) below. 

For s > 3/2, the Neumann condition has to be taken into account, and we 
have in particular 

H%(fl) = {u G H 2 (il) ; u„ = on <9tt} 

which will play a particular role in the sequel. More generally, a similar charac- 
terization holds for 3/2 < s < 7/2. For s > 7/2, additional boundary conditions 
have to be taken into account. 
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The case s = 3/2 is critical (note that w„|an is not well defined in that space) 
and one can show that 

Fj(fi) = LeH?(Q); J^-^-dx <oo\ 
where d{x) denotes the distance to 9fi. A similar result appears in [17]; more 

3 

precisely, such a characterization of H^(Q) can be obtained by considering 
functions u such that u x <E Vo(^) where Vo(fi) is defined in [17] as the equivalent 

1 /2 

of our space (f2) with Dirichlet rather than Neumann boundary conditions. 
We do not dwell on this issue since we will not need this result in our proof. 

3.2 The operator / 

As it is explained in the Introduction, the operator / is related to the computa- 
tion of the pressure as a function of the displacement. From this point of view, 
the operator / is a Dirichlet-to-Neumann operator associated with the Lapla- 
cian. Since we study the problem in a periodic setting we explained that this 
yields to consider Neumann boundary conditions on a cylinder C = fix (0, +oo). 

Spectral definition. It is convenient to begin with the spectral definition of 
the operator /: With and (pk defined by (fT9]) . we define the operator 

oo OO ^ 

I:^2c k tp k i — > -^c fe Af^ fc (20) 

k=Q k=0 

which clearly maps iJ 1 (fi) onto L 2 (Q) and (Q) onto 7J 1 (fi). 

Dirichlet-to-Neuman map. With the spectral definition in hand, we are 
now going to show that / can also be defined as the Dirichlet-to-Ncumann 
operator associated with the Laplace operator supplemented with Neumann 
boundary conditions. 

To be more precise, we consider the following boundary problem in the 
cylinder C = Q x (0, +oo): 

-Aw = in C, 

v(x,0) = u(x) on Q, (21) 
v v = on dfl x (0, oo). 

We will show that we have 

I(u) = d y v(;0). 

We start with the following result which show the existence of a unique harmonic 
extension v: 

i 

Proposition 1 (Existence and uniqueness for (|2T]) ). For all u 6 H^(Q,), there 
exists a unique extension v £ 7J 1 (C) solution of 1(21]). 
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Furthermore, if u(x) = YluLi Ck<Pk(x), then 

oo 

. ^ 1 

v(x,y) = 2_^c k ip k (x)exp(-\£y). (22) 



k=l 



Proof. We recall that H^(Q) = (O), and for a given (Q) we consider 

the following minimization problem: 

inf |y \Ww\ 2 dxdy; we H l (C), w(-,0) = u on o\ . 

Using classical arguments, one can show that this problem has a unique mini- 
mizer v (note that the set of functions on which we minimize the functional is 
not empty). This minimizer is a weak solution of (|21[) . In particular, it satisfies 

Vu • Vio dx dy = 



c 



for all w € i/ 1 (SI) such that w(-,0) = on f2, which includes a weak formulation 
of the Neumann condition. 

Finally, the representation formula (f22|) follows from a straightforward com- 
putation. Indeed, we have 

/ / \Vv\ 2 dxdy = \d x v\ 2 + \d y v\ 2 dx dy 

Jo Jn Jo Jn 



2V c 2A fe / cxp(-2A 

oo 1 
fe=l ^ A fc 

E^ /2 - IL '" 2 



l- /2 y) °^ 



11 . i 

"ff2(!l) 

fe=l 

which shows that w belongs to H l {C). The fact that u satisfies (|2Tj) is easy to 
check. □ 

We can now show: 

Proposition 2 (The operator / is of Dirichlet-to- Neumann type). For all u € 

ifjy(f2), we /lave 

I(u)(x) = — — (x, 0) = d v v(x, 0) /or aZZ a; £ f2, 

where v is the unique harmonic extension solution of (21]) . 
Furthermore I o I{u) = —Au. 
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Proof. This follows from a direct computation using ([212]) . Furthermore, if u is 
in Hff(Cl), we know that J2kLo c k^k < 00 ■ ^ is now easy to derive the following 
equality 

00 

I(I(u)) = ^CkXkipkix) = -Ait. 

fe=0 

□ 



Integral representation. The operator / can also be represented as a sin- 
gular integral operator. Indeed, we will prove the following 

Proposition 3. Consider a smooth function u : f2 — > K. T/ien /or a// x G f2, 

= / (u{y) - u{x))v{x 1 y)dy 
J a 

where v(x, y) is defined as follows: for all x,y G f2, 

y) = 2 (l - cos(7r(a; - j/)) + 1 - cos(^(x + yj) ) ' (23) 

Proof. We use the Dirichlet-to-Neumann definition of /. Let v denote the so- 
lution of (HJ). Then v is the restriction to (0,1) of the unique solution w of 
(|2ip where f2 is replaced with (—1,1) and u is replaced by its even extension to 
(— 1, 1). In particular, w is even with respect to x. Then there exists a holomor- 
phic function W defined in the cylinder ( — 1, 1) x (0, +00) such that w = Re(W). 
Next, we consider the holomorphic function z e lwz = e~' nv e l ' KX defined on the 
cylinder (— 1, 1) x (0, +00) with values into the unit disk D\ = {(x, y) : x 2 + y 2 < 
1}. If 2 denotes the complex variable x + iy, then a new holomorphic function 
Wo is obtained by the following formula 

W(z) = W (e" z ). 

In particular, Wq is defined and harmonic in D\. This implies that the function 
Wq can be represented by the Poisson integral. Precisely, 

i-\z\ 2 r w {Y) 

IdDi 



This implies that for all z e C, 

v ' 27T J_ x \e md - e -*ye lnx \ 2 

and we finally obtain 

1 - e- 2 -v f 1 w(0,O) 

Wlx.y) = / 1 - a : 7^ dU. 

v >»J 2 J_ 1 \e™ e - e -^Ve mx \ 2 
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Taking w = 1, we get in particular the following equality: 

i - e- 2 *y f 1 1 
1 = / — — a — -de. 

We deduce: 

w(x,y)-w(x,0) l-e~ 27ry f 1 w(0, 0) - w(x, 0) 



d0 



which implies (letting y go to zero): 

d y w(x,o) = tt J ^ -^—^——^de. 

The integral on the right hand side of the previous equality is understood in the 
sense of the principal value of the associated distribution. We finally use the 
fact that w is even in x and is equal to u on D, to obtain the following singular 
integral representation of I(u): 

1 1 1 



□ 



I(u)(x) = it I {u{6, 0) - u(x, 0)) ( |1 _ eiT(a _ 9)|2 + |1 _ eWa+e)|2 ) dB. 



The space H-?(n). The space H ? (f2) is defined as the topological dual 
space of H ? (57). It is classical that for any u E 2 (£7) ; there exists ui £ L 2 (f2) 
and it 2 £ H^{Q) such that u = u\ + d x U2 (in the sense of distributions). We 
will also use repeatedly the following elementary lemma, whose proof is left to 
the reader: 

Lemma 1. If u £ Hi(Q), then the distribution I(u) is in H~i(£l) and for all 

v £ H^(n), 

+00 ^ 

(I(u),v) 1 1 =— > X? Chdh 

k=0 

where u = J^t^o c Wk and v = J^t^o dk^k- In particular, 
— (I(u),u) 1 1, = Hull 2 1 

Important equalities. The semi-norms || • \\^ (Q) , II • lljji(n), II ' ll^f (n) 
and || • 1 1 jj 2 ,q% are related to the operator / by equalities which will be used 
repeatedly. 

Proposition 4 (The operator I and several semi-norms). 
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For all u G H 2 (fi), we ftawe 



1 



9 JnJn^ ~ U<yV ^ v ( x,y } dxdy = 
For m G H 1 ^), we have 

(T(u))' 2 dx = \\u\\% 1(n) 
For all u G H^(fl), we have 

I(u) x u x dx = \\u\\ 2 3 
For all u G H^(fl), we have 

r (d x I(u)) 2 dx=\\u\\% lm . 
Remark 1. Note that I(u) x ^ I{u x ). 

Proof. The two first equalities are easily derived form the definition of /, def- 
initions of the semi-norms, the integral representation of / and the fact that 
v(x,y) = v{y,x). 

In order to prove the third and fourth equalities, we first remark that d x Lp k = 
1 

— A| sin(fc7rx) form an orthogonal basis of L 2 (Q). 

In order to prove the fourth equality, we first write 



d x (I(u)) = - CfcAf d x cpk inL 2 (0) 

from which we deduce 

/ (I(u) x ) 2 dx = Vc^Afe / (8 x ^ k ) 2 dx 



fc=i 



00 „ 

y^c^Afc / tp k (-d xx (p k )dx 
fe =i Jn 

OO 

^2 c l x l = \ \ u \\ 2 hi- 



fc=0 

As far as the third equality is concerned, we note that 



u x = 

k=0 
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We then have 



c 00 1 c 

- \ I{u) x u x dx = ^c|A| I (d x (fi k ) 2 dx 

«/ O q J a2 

00 if 

fc=o ^ a 

= Xl c fe A fc / ^wldx 
k=0 Jn 

00 

= Y4\l =\\u\\ 2 3 . 



fc=0 



□ 



3.3 The problem —I(u) = g 

We conclude this preliminary section by giving a few results about the following 
problem: 

For a given g g £ 2 (fi), find it 6 i/ 1 (fi) such that 
-I(u)=g. 

Note that J n J(u) da; = for all u 6 (since J n ifikdx = for all > 1) 

and so a necessary condition for the existence of a solution to ([24)) is 

(7(2:) da; = 0. 

a 

Note also that there is no uniqueness since if u is a solution then u + C is also a 
solution for any constant C. We may however expect a unique solution if we add 
the further constraint / u dx = 0. Indeed, a weak solution u £ H~% (fi) for g € 
H~i(Vl) can be found using Lax-Milgram theorem in {u € Hi (ft) ; J n udx = 
0} equipped with the norm ||«|| . 1 . Alternatively, we can use the spectral 

framework: For g <G L 2 (Q) such that J Q g{x) dx = 0, we have 

00 00 
g = dk^k with d\ < 00. 

k=l k=l 

We can then write: 

00 , 

u = r\g):=J2^k (25) 
*=i ^ 

which clearly lies in iJ 1 (J7) and satisfies J n udx = 0. The fact that the ipk's 
form an orthogonal basis of L 2 (Q) implies that there is only one solution of (|24p 
such that J Q udx = 0. Finally it is clear from (f25j) that further regularity on g 
will imply further regularity on u. We sum up this discussion in the following 
statement 
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Theorem 3. For all g G L 2 (fl) such that J^gdx = 0, there exists a unique 
function u G H 1 ^) such that —I(u) = g in L 2 (Q) and J^udx = 0. Further- 
more, if g is in H 1 ^), then u G H^(Q). 

We will also use the following corollary of the previous theorem 

Corollary 1. For all g G L 2 (Q), there exists a unique solution u G H l (£l) of 

—I(v) + / v dx = g. 
Jn 

Proof. Set m = J n g(x)dx and consider g = g — m. Then g € L 2 (f2) and 
J n gdx — 0. There is a (unique) u G iJ 1 (SI) such that 

—I(u) = g — m, / u(x) dx = 0. 
Jn 

We then set v = u + m. Then J^v dx = m and 

—I(v) = —I(u) — g — m = g — / v dx. 



n 

As far as uniqueness is concerned, if we consider two solutions v\ and vi then 
we have 

/ v\dx = / V2dx = g 
Jn Jn Jn 

and this implies that w = v\ — v-i satisfies —I(w) = 0. The uniqueness of the 
solution given by Theorem [3] implies that w = and the proof is complete. □ 

4 A regularized problem 

We now turn to the proof of Theorem [T] The degeneracy of the diffusion co- 
efficient is a major obstacle to the development of a variational argument. As 
in the existence of solution for ([5]) is thus obtained via a regularization 
approach: Given e > 0, we consider 

dtu + d x (f £ (u)d x I(uj) = 0, te{0,T),xen (26) 

where 

f £ {s)= S+ n +s 
(with s+ = max(s,0)), with the initial condition 

u(0,x) = u {x) (27) 

and boundary conditions 

u x =0, f £ (u)d x (I(u))=0 ondn. 
The first step in the proof of TheoremQ] is to prove the following proposition: 
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Proposition 5 (Existence of solution for the regularized problem). For all 

uq G H^(n) and for all T > 0, there exists a unique function u £ such that 

u £ G L°°(0, T; Hi (fi)) n £ 2 (0, T; f$(fi)) 

solution of 

u e d t tpdxdt+ // f e (u £ )d x I(u £ )d x ipdxdt = — / uoVj(0, •) da; (28) 
i Jjq Jq 

for all tp G ^([0,^,^(0)) w't/i Q = x (0,T). 
Moreover, the function u £ satisfies 

u £ (x,t)dx= / w (a;)da; o.e. fe(0,T) (29) 



t 

e\\2 



and 

H ue (*'-)llkro. +2 / / /«(« e )(0*/(« e )) a <fcds< Nf. i . o.e. te(0,T). 

(30) 

Finally, if G e is a non-negative function such that G"(s) = j ^ , then u £ 
satisfies for almost every t G (0, T) 



G e (u £ ){x,t)dx+ / |U s (s)||a ds< / G £ (u )dx. (31) 
Jo ^ (n) Jn 

Remark 2. Note that this result does not require condition (J8j) to be satisfied 
and is thus valid with compactly supported initial data. However, we will need 
condition (J5J to get enough compactness on u £ to pass to the limit e — > and 
complete the proof of Theorem [TJ 

There are several possible approaches to prove Proposition [5] In the next 
sections, we present a proof based on a time discretization of (|28p and fairly 
classical monotonicity method (though the operator here is not monotone, but 
only pseudo- monotone). 

4.1 Stationary problem 

In order to prove Proposition [5j we first consider the following stationary prob- 
lem (for r > 0): 

For a given g G H 2 (O), find u G H^(Jl) such that 
f u + rSxC/eCuja^JCu)) = .g in fi (32) 
I d x u = and d x I(u) = on 



Once we prove the existence of a solution for (|32[) . a simple time discretization 
method will provide the existence of a solution to (|28|) . We are going to prove: 
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Proposition 6 (The stationary problem). For all g e H^(Q), there exists 
u e Hff(Q) such that for all ip e i7 1 (fi), 

- I {u — g)ipdx- [ f e (u)d x I(u)d x ipdx = 0. (33) 



Furthermore, 

u(x) dx = I g(x) dx , (34) 



\W\\% hm +^JJ e {u){dJuYdx < \\gf khn) , (35) 

and if G e (g) dx < oo then 

{ G e {u)dx + r\\u\\ 2 3 < [ G £ (g)dx. (36) 
A2 h$(ci) Jn 

In order to prove such a result, we have to reformulate ([33]): 

New formulation of (|33j) . We are going to use classical variational methods 
to show the existence of a solution to (J33J. In order to work with a coercive 
non-linear operator, we need to take ip = —I(v) as a test function. We note, 
however, that by doing that we would restrict ourself to test functions with 
zero mean value. In order to recover all test functions from H (Cl), we use 
Corollary [T] and consider 

ip = -I(v) + vdx (37) 



n 



for some function v £ H]^(Cl). Let us emphasize the fact that Corollary[T]implies 
in particular that there is a one-to-one correspondence between ip G i/ 1 (J7) and 
v G H%(Q) satisfying ([37] , 

Using (|37)) , Equation ([33]) becomes: 



n 



ul(v) dx + ( / udx] ( / v dx j +t j f e (u)d x I(u)d x I(v) dx 

g I(v) dx + (J g dx^j (J v dx^j . (38) 

We can now introduce the non-linear operator we are going to work with. 

A non-linear operator. We define for all u and v G i/jy(f2) 

A(u)(v) — — / ul(v)dx+i / udx) I / vdx)+r f e (u)d x I{u)d x I(v)dx. 
Jn \Jn J \Jn J Jn 

One can now show that this non-linear operator is continuous, coercive and 
pseudo-monotone. Classical theorems imply the existence of a solution to the 
equation A(u) = g for proper g's. More precisely, we have the following propo- 
sition: 
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Proposition 7 (Existence for the new problem). For all g £ (f2) there exists 
u G Hfj(fl) such that 

A(u)(v)=- gI(v)dx+( gdx) ( I vdx) for all v G H%(fl). (39) 
Jn \Jn J \Jn ) 

For the sake of readability, we postpone the proof of this rather technical 
proposition to Appendix and we turn to the proof of Proposition [6] 

Proof of Proposition^ For a given g G H^(fi), Proposition [7] gives the exis- 
tence of a solution u G Hff(fl) of (|38|) . We recall that for any cp G H 1 ^), there 
exists v G Hfj(rt) such that 

<p = —I(v) + / vdx 
Jn 

and so equivalently, we have that u satisfies (f3"3"]) for all ip G H 1 ^). 

Next, we note that the mass conservation equality p4p is readily obtained 
by taking v = 1 as a test function in (|3"5)l , while ([33)) follows by taking w = 
u — J n udx: 

h\\ 2 . hm +r [ f e (u)\d x I(u)\ 2 =- [ gl(u)dx 
Jn Jn 

^llffll^ ( n)NI^(n)^2 llff|l ^(n) + 2 IMI ^cn)- 

Finally since G' £ is smooth with G' £ and G" bounded and fi is bounded, we 
have G' e (u) G H l (£l). We can thus find v G H%(Q) such that 



n 



uG' e (u) dx + r f e (u) F'J{u) d x I{u) d x u dx = - gG' e (u) dx 
n Jn Jn 

and so (using the definition of G E given in Proposition [5]) 

— t / d x I(u) d x udx = / G' £ (u)(g — u) dx 



-I(v)+ I v(x)dx = G' s (u). 
Equation (|38|) then implies: 



Since G e is convex (G" > 0), we have G' e (u)(g — u) < G e (g) — G e (u) and we 
deduce (|36|) (using Proposition [4}. □ 
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4.2 Proof of Proposition [5] 

In order to construct the solution u £ of (f2"tj|) , we discretize the problem with 
respect to t, and construct a piecewise constant function 

u T (x,t) =u n (x) forte {tit, (n + l)r),ne {0, . . . , N + 1}, 

where r = T/JV and (w n )ne{o,...,./v+i} is such that 

-{U n+1 - U n ) + d x {fe{u n+1 )d x I{u n+1 )) =0. 
T 

The existence of the u n follows from Proposition [5] by induction on n. We 
deduce: 

Corollary 2 (Discrete in time approximate solution). For any N > and 
Uq € if5(f2), £/iere exists a function u T £ L°°(0, T; _ff 3 (17)) suc/i £/ia£ 

• 1 1 — ^ m t (x, i) is constant on [fcr, (fc + l)r) /or fc S {0, . . . , iV}, T = j, , 

• w r = uo on [0, r) x fi, 

• for all ip e C^O,!^ 1 ^)), 

'J—-^-tpdxdt= [[ f E (u T )d x I(u T )d x <pdxdt (40) 

where Q t ,t = (t, T) x 17 and S T u T (x,t) — u T {t — t,x). 
Moreover, the function u T satisfies 



u T {x,t)dx= I u (x)dx a.e.te{0,T) (41) 

! in 

and /or aZ/ 1 £ (0, T) 

IKfcOllli, +2 f I fe{u T ){d x I{u T )) 2 dxdt < \\uo\\ 2 i (42) 
and z/ J n G £ (mo) dx < oo, then for all t € (0,T) 

/ G 6 (u r (i,-))dz + / ||u T || 2 3 ds < / G e (u )cfo. (43) 
Jo h%(u) Jq 

It remains to prove that u T converges to a solution of (|28|) as r goes to 
zero. This is fairly classical and we detail the proof for the interested reader in 
Appendix [SJ 
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5 Proof of Theorem [D Case n > 2 



Proposition [5] provides the existence of a solution u £ 6 L°°(0, T; H 2 (Vt)) n 
L 2 (0,T;H%(fl)) of ([28]). Our goal is now to pass to the limit e ->• 0. We 
point out that at this point, the solution u e may change sign and that it is only 
at the limit e — > that we are able to recover a non-negative solution, using the 
fact that n > 2. 



Step 1: Compactness result. First, we note that ([50)1 implies 

IKWiy (n) < IK(*)iy (n) forall £ >0. (44) 

The bound (|44|) and Sobolcv embedding theorems imply that the sequence 
(u £ ) e> o is bounded in L°°(0, T; L p (f2)) for all p < 00 and so f e (u £ ) is bounded in 
L°°(0,T;LP(n)) for all p < 00. Furthermore, j3DJ also gives that f e {u £ )^d x I(u e ) 
is bounded in L 2 (0, T; L 2 (0)). We deduce that 

f e (u £ )d x I{u e ) is bounded in L 2 (0, T; L r (fl)) 

for all r £ [1,2). Writing 

d t u £ = -d x {f e {u s )d x I{u s )) in 2?'(0), 

we deduce that {d t u e ) 6>0 is bounded in L 2 (0, T; W' 1 ^ (CI)) for all r' G (2, +00). 
Thanks to the following embeddings 

#4(0) i 9 (0) -> W" 1 ' r '(0) 

for all q < 00, if follows (using Aubin's lemma) that (u 6 ) e >o is relatively compact 
in C°(0, T, L 9 (il)) for all q < +00. Hence, we can extract a subsequence, still 
denoted by u e such that 

u £ — > u in C°(0, T, £«(fi)) for all g < 00 

and 

u £ — > u almost everywhere in Q. 



Step 2: Passing to the limit in Equation (|28[) . We now have to pass to 
the limit in ([25]). We fix tp G U(Q). Since u e -> it in C°(0, T, L^fi)), we have 



u £ dttpdxdt^ II udftpdxdt. 
Next, we remark that (|30|) implies 



(d x I(u £ )) 2 dxdt<±\\u \\ H i {ny 
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Cauchy-Schwarz inequality thus implies 

ed x I(u e )d x (p dx dt 



Finally, ^ implies that (u £ )ld x I{u £ ) is bounded in L 2 (0, T, L 2 {£1)). Since 
u e is bounded in L°°(0, T; L p (f2)) for all p < oo we deduce that (u £ )™ d x I{u e ) is 
bounded in £ 2 (0, T; £ r (S!)) for all r <G [1,2) and so 

h e := (u e )ld x I(u E ) in L 2 (0, T; Z/(fi))-weak. 

Passing to the limit in (|28j> . we get: 



// udt<pdxdt+ // hd x (pdxdt = — // UQip(0,-)dxdt 
JJq JJq JJq 

for all p € X>(Q). 

Step 3: Equation for the flux h. In order to get ||7J), it only remains to 
show that 

ft = u n + d x I{u) 

in the following sense: 

h<j)dxdt = - nu1~ 1 d x u I(u) <j)dxdt - u\l (u) d x (j) dx dt (45) 
q JJq JJq 

for all test function <f> such that (j>\dQ. = 0; that is 

h = d x {u\I{u)) - nv%T l {d m u)I{u) in V'(Cl). 
For that we note that since 

G e (uo) dx < C, 

Inequality ([3Tj) implies that iu e ) £ is bounded in L 2 (0, T; Hi (fi)). Recall that 
{d t u € ) e is bounded in L 2 (0, T, M / ~ 1,r '(r2)) for all r' € (2,+oo). Aubin's lemma 
then implies that u e is relatively compact in L 2 (0,T; H s (ft)) for s < 3/2. In 
particular, we can assume that 

I(u e ) — > I(u) in L 2 (0,T;L 2 (fi)) 



and 



d x u £ — ► d x u in L 2 (0,T;L p (fi)), for all p < oo. 



Writing 

/i e </> = I I (u e )ld x I(u £ ) <f)dxdt 



Q 



n(u e )T'd x u e I{u e ) <j)dxdt- / / (u e )+ d x <f>dxdt 



we see that those estimates, together with the fact that u e converges to u in 
L°°(0, T; L p (fl)) for all p < oo, are enough to pass to the limit and get (05]). 



24 



Step 4: Properties of u. It is readily seen that u satisfies the conservation 
of mass (fT0|) (by passing to the limit in (j29|0 . and the lower semicontinuity of 
the norm implies the entropy inequality (|12[) . 

Next, Inequality implies that g £ = {u £ + )~z d x I{u £ ) converges weakly in 
L 2 ((0, T) X f2) to a function g, and the lower semicontinuity of the norm implies 
(fTT]) . Proceeding as above we can easily show that 

g = d x {ull{u)) - ^u\- l d x u I(u) in £>'(Q). 

Step 5: non-negative solutions. It remains to prove that u is non-negative. 
This will be a consequence of (f5T|) and the fact that n > 2. Indeed, we recall 
that for all t we have 

f G E (u £ (t)dx< [ G e {u )dx (46) 
Jn Jn 

where is such that G''(s) = -, — }„ , . As noted in the introduction, we can take 



G e (s) = G'l(t)dtdr 



which is a nonnegative convex function for all e. Noticing that we can also write 
G e (s) = f f G"(t) dtdr when s < 1, it is readily seen that G E (s) is decreasing 
with respect to e (so G e {s) < Gq(s) for all e > 0). Hence 



We deduce (using (|46 



G e (uo)dx< / Gq(uo) dx < +oo . 



limsup / G e (u e {t))dx < +oo . (47) 
e-^o Jo 



Next, we recall that u e (-,t) converges strongly in L p {Vl) to u(-,t). We can 
thus assume that it also converges almost everywhere. Egorov's theorem then 
implies the existence of a set A v C H such that u E (-,t) — > u(-,t) uniformly in 
A n and |Q\ A n \ < i]. For some 5 > 0, we now consider 

C n ,s = A, n {«(•,*) < -2d}. 

For every 77, S > there exists £0(77, <5) such that if e < eo(r],5), then u £ (-,t) < 
-8 in C, )5 . 

But this implies that C Vt s has measure zero. Indeed, if not, then for e < 
£0(77, 6) we have 

G e (u £ (x,t)) > G £ {-6) — > G (-5) = +00 for all x G C n j 

(we use here the assumption n > 2) and by Fatou lemma, we get 

liminf / G e (u e (x, t)) dx > / liminf G e (u £ (x, t)) dx = +00 

e->o I n £->0 



25 



which contradicts (|Tf)) . 

We deduce that for all 8 > and all 77 > we have 

< -2«5}| < ia ; , 5 | + 1«\ a,i < ?? 

and so \{u(-, t) < — 2S}\ = for all S > 0. We can conclude that 



K,t)<o}= |J /«(.,*)< -I] 

n>l ^ 



has measure zero and so u[x,t) > a.e. x G and for all i > 0. 

6 Proof of Theorem [D Case n 6 [1,2) 

When n G [1,2) the entropy inequality cannot be used to prove that lim £ _>o u e 
is non-negative. We thus proceed as in Bertozzi and Pugh |14j : Introducing 

fs(s) 



5s n + s 3 
and 

Uq(x) = u (x) + S. 

For n < 2, we have fs(s) ~ s 3 /<5 as s — > 0, and so the corresponding entropy 
Gg, defined by 

satisfies 

Gs{s) = d(± + ^-l)+G {s) 

where Gq(s) is bounded in the neighborhood of s = 0. It is thus readily seen 
that there exists C such that 

G e (u s (x))dx < C. 

n 

Furthermore, we have Gs(0) = +oo, so the proof developed in the previous 
section (regularizing the equation by introducing fs.e(s) — fs(s)+s) implies the 
existence of a non-negative solution u of 

d t u s + d x (f s (u s )d x l(u 5 )) = 

satisfying the usual inequalities (mass conservation, energy and entropy inequal- 

ity). 

Proceeding as in the previous section, we can now show that the sequence 
(u s )s>o is relatively compact in C°(0, T, L q (Cl)) for all q < +oo and in L 2 (0, T; H s (fl)) 
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for s < 3/2. In particular, we can extract a subsequence, still denoted u s , such 
that 



U S — > u 



in C°(0,T,L«(Q.)) for all q < +00 
u s — > u almost everywhere in Q 
I(u s )-^I{u) in L 2 (0,T;L 2 (Q)) 
d x u s — > d x u in L 2 (0, T; L p (n)), for all p < +00. 

Furthermore, since u s > for all 5 > 0, we have 

u > a.e. (a:,i) € x (0,T). 



In order to pass to the limit in the equation, we mainly need to check that 
fs(u s ) (respectively /^(u 5 )) converges to u n (respectively nu 71 ^ 1 ) in L P (Q). This 
is a direct consequence of the convergence almost everywhere of u s , Lebesgue 
dominated convergence theorem and the fact that 

f 5 (s) < s n for all s > 

and 

o r 2+2n , „5+n 

f's(s) = (£5 , t + + s " 2 < ^ + 2)s^ for all s > 

(note that we need n > 1 to complete this computation). 
This complete the proof of Theorem [TJ 



7 Proof of Theorem [2] 

In this section, we prove Theorem [2j We recall that n > 3 and we consider the 
sequence u e of solution of the regularized equation (|26l) introduced in the proof 
of Theorem [TJ 

We recall that inequality (f3"Tj) implies that (ir) e is bounded in L 2 (0, T; (O)). 
Since (5 t w e ) e is bounded in L 2 {0,T,W~ iy (n)) for all r' G (2, +00), Aubin's 
lemma implies that u e converges strongly in L 2 (0, T; C a (J7)) for all a < 1. We 
can thus find a subsequence such that u e (t) converges strongly in C a (£l) for 
almost every t (that is for all t G P, where |(0, T) \ P| = 0). 

Next, we note that for t G P, it is actually strictly positive. Indeed, if 
u(xojio) = 0, then for any a < 1, there is a constant C a such that 

u(x,t) <C\x-x \ a 

We deduce 

j G{u{x,t Q ))dx> J {Ca]x _ 1 Xo|Q)n _ 2 ^ 
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Given n > 3 we can choose a < 1 such that a(n — 2) > 1. We deduce 

G(u(x, to)) dx = oo 



which contradicts ([3"TT) . 

We deduce that there exists 6 > (depending on i) such that for e small 
enough 

u £ (-,t) > 6 in Q. 

Next, we note that after removing another set of measure zero to P. we can 
always assume that 

liminf / Uu s )\d x I(u E )\ 2 dx < oo for all t G P. 

E->0 J 

Indeed, if we denote 

A k = {t G P ; liminf / / £ (u e )|a 2 J(u e )| 2 rfa; > k} 
we have (using Fatou's lemma): 

C > liminf ^ J f e (u e )\d x I{u e )\ 2 dx dt 



> liminf / / / e (u e )|0*J(ti e )l 2 dar<ft 



> / liminf / f e (u e )\d x I(u e )\ 2 dxdt 
JA k J 

> k\A k \ 

We deduce \A k \ < C/k and thus 

\{teP; liminf / /efu^l^/fu 6 )! 2 = oo}| = 0. 

It follows that for t G P, we have 



and so 



liminf |9 2 -/(u e )| dx < oo 



(•,i) u(-,t) in ifjy- (fi)-weak. 



In particular we can pass to the limit in the flux J e = f e (u £ )d x I(u e ) and 
write 

lim ,L = J= f(u)d x I(u) in L^Q), a.e. t G (0,T) . 

£->0 

Furthermore, we note that we recover the boundary condition in the classical 
sense: 

u x (x,t) = for x G dtt and a.e. £ G (0, T). 
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A Proof of Proposition [7] 

We denote 

V = H%(Q.). 

For any u € V the functional A{u) is clearly linear on V and since V is contin- 
uously embedded in L°°(il), we have 

\A(u)(v)\ < [lMI H4(n) +r(e + ||tt||^)||«||v] |H|v. (48) 

(Note that we used Proposition@]to get this inequality) . The non- linear operator 
A is thus well-defined as a map from V to V. Moreover, it is bounded. 
Next, we remark that we have 

2 

2 



A(u)(u)>- / uI(u)dx+[ / udx) +e \d x I(u)\ 2 dx. 
Jn \Jn J Jn 

We deduce from Proposition |4] that 

A{u){u)> T e\\u\ 2 Hl(n) . (49) 

In particular, we have 

A{u){u) 



" ■ v 



— >• +oo as M y — ^ +oo. 



The operator A is thus coercive. Proposition [7] will now be a consequence of 
classical theorems if we prove that A is a pseudo-monotone operator. Since we 
already know that A is bounded, it remains to prove the following lemma: 

Lemma 2 (A is pseudo-monotone). Let Uj be a sequence of functions in V such 
that uj u weakly in V . Then 

liminf A(uj)(uj — v) > A(u)(u — v). 

j 

Before we prove this lemma, let us notice that for g € iJs (57) ; the application 

T g : v - / gl(v) dx + ( / g dx I ( / v dx I 

belongs to V'. Hence, using Theorem 2.7 (page 180) of [31], we deduce that for 
all g G (fi), there exists a function tt € V such that A(u) = T g in V', which 
completes the proof of Proposition [7] 



It remains to prove Lemma [5J 
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Proof of Lemma\^ We first write 



A(uj)(uj — v) = Ujl(uj — v) dx + I / uj dx I I / (uj—v)dx 

Jn \Jq J \Jn 

+ T / f e (uj)d x I(Uj)d x I(uj-v)dx 
Jn 

= II Uj II 2 i — (in, v) i 
11 JU H?(n) v J ' 'hi 



where 



+r / /eCujO^AujOr-T / f e ( Uj )d x (I Uj )d x (Iv)(50) 
Jn Jn 



(u,v) i =— / ul(v) dx + I / uefec) I / u dx 
Jn \Jn J \Jn 



We need to check that we can pass to the limit in each of those terms. 
Since Uj converges weakly in V we immediately get 



and 



liminf ||u,-|| i ■ ||w|| 2 i 
j-v+oo JU H2(n) - "H2(n) 



j ]hnJu j ,v) Hi =-{u,v) Hi 



Since Uj is bounded in £fjy-(fi), it is compact in L°°(f2), and so f e (uj) converges 
to f e (u) strongly in L 00 (d). We thus write 

f e ( Uj )(d x i( Uj )f = f (f £ (u 3 ) - f £ {u)){d x i{u )f + f f £ (u)(d x i( Uj )) 2 

n Jn Jn 

> -\\fe(u j )-fe(u)\\ L - ( n ) \\u j \\ 2 v + [ f e {u){d x I{u 3 )f. 



n 



The first term goes to zero and we have 

V r f £ Jujd x I(u 3 ) -± V f f £ Jujd x I(u) in L 2 {tt). 
Again, the lower semicontinuity of the L 2 -norm gives 

lim t f f £ (u)(d x I(u 3 )) 2 > f f £ (u)(d x I(u)) 2 . 
3^°° Jn Jn 

Finally, we have 

f £ (uj)d x I(v) -> f £ (u)d x I(v) in L 2 (f2) strong, 
d x I(uj) — ^ d x I{u) in L 2 (i7) weak 

which gives the convergence of the last term in (f5U|) and completes the proof of 
the lemma. □ 
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B Proof of Proposition [5] 

The proof of Proposition [5] is divided in three steps. 



Step 1: a priori estimates. We summarize the a priori estimates in the 
following lemma: 

Lemma 3 (A priori estimates). The solution u T constructed in Corollary \^ 
satisfies 



I i < Kll i , (51) 

>/S||0x/(u T )IU»W) < C, (52) 
- S T u T 

T 

for all r' G (2, +oo) where C does not depend on t > (but does depend on r') 



< C, (53) 

i 2 (r 1 T 1 W- 1 .' r '(n)) 



Proof. Estimate (jolj) and (|52|) arc direct consequences of (|41|) and (|42 
Next, we note that 

The bound (f5"Tj) and Sobolev embedding theorems imply that the sequence 
(ii T ) r>0 is bounded in L°°(0, T, L p (f2)) for all p < oo and so f e (u T ) is bounded 
in L°°(0,T,LP(f2)) for ah> < oo. Since d x I{u T ) is bounded in L 2 (Q), we deduce 
that f e {u T )d x I{u T ) is bounded in L 2 {t,T, L r {fl)) for all r G [1,2). It follows 
that 

d x {f e {u T )d x I{u T )) is bounded in L 2 (t,T, W^' '(fi)) 
for all r' G (2, oo). □ 

Step 2: Compactness result. Thanks to the following imbeddings 

F5(0) <^> L 9 (Q) -> ^-^'(Q) 

for all q < oo, we can use Aubin's lemma to obtain that (u T ) T is relatively 
compact in C°(0, T, L 9 (fi)) for all g < oo. 

Remark that {d x I{u T )) T is bounded in L 2 (Q) and (m t ) t is bounded in 
L^O.TjL 1 ^)). It follows that (w T ) r is bounded in L 2 {Q,T,H 2 N {Vt)). Since 

3 

we deduce that {u T ) T is relatively compact in L (0,T; H^(fl)). Up to a subse- 
quence, we can thus assume that as r — > 0, we have the following convergences: 

• u T — y u e G L°°(0, T, H =(0)) almost everywhere in Q; 

• u T ^u e in L 2 (0, 7,^(0)) strong; 

• d x I{u T ) — ^ d x I(u e ) in L 2 (Q) weak. 
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Step 3: Derivation of Equation (|2"8|) . We want to pass to the limit in (|4()p. 
We fix g £*([(), T), if 1 (Q)). Then 

u T - S T u T ff ip(x,t) - ip(t + T,x) 
<P= // « OM) 



— / / u T (x, t)(p(x, t) dx H — / / u T (x,t)ip(x,t) dx. 
T Jo Jq t Jt-t Jq 



We deduce: 

u T — S T u T 
T 



<p->- !! u £ (d t <p)- I u £ (0,x)ip(0,x) 
J Jo Jq 



dx + 0. 



It remains to pass to the limit in the non-linear term. Let r\ > 0. Since u T — > u e 
almost everywhere in Q, Egorov's theorem yields the existence of a set A v C Q 
such that \Q \ A n \ < ?y and 

u T — > u £ uniformly in . 

In particular, 

Vfe(u T )d x <p -> y/f s (u e )d x <p in L 2 (A V ) 

and 

y/M^)d x I(u T ) - ^Jj^)d x I{u £ ) in L 2 {A V ) . (54) 

Hence 



f e {u T )d x I{u r )d x ip -> / f £ (u e )d x I{u £ )d x (p 
as r goes to zero. 

Finally, we look at what happens on Q \ A v . Choose Pi,P2,P3 such that 
J^p" 1 = 1 and write 



\f e (u T )d x I(u T )d x ip\ 

Q\A V 

< ||5 x ^|| L =o (0 . T . LP1(n)) / \\fe(u T )d x I{u T )\\ LP2in) \\l Q \ A J\ LP3in) 
Jo 

< ll^^llioofo ,T,LPi(n)) II fe(u T )d x 7(w T )|| L 2( ,T,LP2(Q)) II lq\A„ II L 2 (0,T,LP3(n)) 

We now choose P2 € [1, 2) (and so pi > 2 and pz > 2). 

|/ c (u T )0 x J(u r )0 x ¥>l < C(^)||1q VA J| LP 3 ( q) < Cfa)^ . 
Since rj is arbitrary, the proof is complete. 
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Step 4: Inequalities. Since u T — > u £ in L°°(0, T, mass conservation 

equation ([29|) follows from (|4T|). 

Estimate (f3"u]l follows from (|4"2"j) . Indeed, since it T — > u e almost everywhere, 
Proposition |4] and Fatou's lemma imply that for almost t <G (0,T) 

\\u e {t)\\ 2 . 1 <liminf ||u T (t)|| 2 . x . 
Thanks to (|54"1) , we also have 

I(u £ )\ 2 dx dt < limM [ I f E (u T )(d x I(u T )) 2 dxdt 

f e (u £ )\d x I{u £ )\ 2 dxdt. 



o Jo Jo Jo 



JO\A„ 

Letting i] — > permits to conclude. 

To derive (pTTj) we note that G E (u T ) — ?> F e (u £ ) almost everywhere. So Fatou's 
Lemma implies for almost every t € (0,T) 

G e (w £ (a;,0)^ < liminf / G ^{u T {x, i)) dx < / G e (u )da;. 
T -* Jn Jn 

3 

Finally, since (u T ) T is relatively compact in L 2 (0, T; H^(ft)), we have 



||u 8 (s)|K 3 ds = lim / s 



and so (I3T1) follows from (l43l). 
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